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ABSTRACT 

The Radiation Assessment Detector (RAD), on board Mars Science Eaboratory’s 
(MSE) rover Curiosity, measures the energy spectra of both energetic charged and 
neutral particles along with the radiation dose rate at the surface of Mars. With these 
first-ever measurements on the Martian surface, RAD observed several effects influ¬ 
encing the galactic cosmic ray (GCR) induced surface radiation dose concurrently: [a] 
short-term diurnal variations of the Martian atmospheric pressure caused by daily ther¬ 
mal tides, [b] long-term seasonal pressure changes in the Martian atmosphere, and [c] 
the modulation of the primary GCR flux by the heliospheric magnetic field, which cor¬ 
relates with long-term solar activity and the rotation of the Sun. The RAD surface dose 
measurements, along with the surface pressure data and the solar modulation factor, 
are analysed and fitted to empirical models which quantitatively demonstrate how the 
long-term influences ([b] and [c]) are related to the measured dose rates. Correspond¬ 
ingly we can estimate dose rate and dose equivalents under different solar modulations 
and different atmospheric conditions, thus allowing empirical predictions of the Mar¬ 
tian surface radiation environment. 


'institute of Experimental and Applied Physics, Christian-Albrechts-University, Kiel, Germany 

guo@physik.uni-kiel.de 

^Southwest Research Institute, Earth, Oceans & Space Department, Durham, NH, USA 
^Southwest Research Institute, Space Science and Engineering Division, Boulder, USA 
"'NASA Headquarters, Science Mission Directorate, Washington DC, USA 
^Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA 
^Einnish Meteorological Institute, Helsinki, Einland 

^ Aerospace Medicine, Deutsches Zentrum fiir Luft- und Raumfahrt, Koln, Germany 



- 2 - 
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1. Introduction and Motivation 


The assessment of the radiation environment is fundamental for planning future human mis¬ 
sions to Mars and evaluating the impact of radiation on the preservation of organic bio- signatures. 
Contributions to the radiat i on environment on the Martian surface are very complex (e.g. JSaganti et al. 
2002 : Dartnell et al. 2007 : Ehresmann et al. 2011 : Kohler et al. 2014 : Ehresmann et al. 2014 ): en¬ 
ergetic particles entering the Martian atmosphere either pass through without any interactions with 
the ambient atomic nuclei, or undergo inelastic interactions with the atmospheric nuclei creating 
secondary particles (via spallation and fragmentation processes), which may further interact while 
propagating through the atmosphere. Einally all primary and secondary particles reaching the sur¬ 
face may also interact with the regolith and, amo ngst others, produce neutrons which could be 
backscattered and detected as albedo neutrons (e.g.. lBovnton et al.ll2004l) . Therefore, the radiation 
environment measured at the surface of the planet is determined by the characteristics of the pri¬ 
mary radiation incident at the top of the atmosphere, the composition and mass of the atmosphere, 
and the composition and density of the surface soil. The above process can be described by a 
simplified mathematical equation: 


P Az, t, E) 


CO 


P-’ Po)Poi(.Eor t)dEQ, 


( 1 ) 


where Eq is the energy of a primary particle with species i (e.g., protons, alpha particles and 
heavy ions) reaching the top of the Martian atmosphere; EqAEq, t) is the spectrum (in the unit of 
counts/MeV/sec/cm^) of primary particle type i at time t; M,y(z, E, Eq) is the yield matrix, repre¬ 
senting the interaction between particles and the atmosphere (and the regolith), of particle type i 
with energy Eq generating particle type j with energy E\ Mij therefore depends on the altitude z 
(or atmospheric pressure E); finally Ej{z, t, E) (in the unit of counts/MeV/sec/cm^) is the result¬ 
ing particle spectra of type j at time t. Compared to Neutron Monitors on Earth which measure 
the count rates of secon dary particles generated by primary fluxes going through the atmosphere 
(Idem & Dormanll2000h . RAD measures a mix of primary and secondary particles. Eurther, there 
is no need to include the geomagnetic cutoff energy in the case of Mars due to the absence of a 
global magnetic field. 
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There are predominantly two types of primary particles (Foi(Eo, t)) reaching Mars: galactic 
cosmic rays (OCRs) and solar energetic particles (SEPs). SEPs are sporadic and impulsive events 
and take place much more frequently during the active phase of the solar cycle. SEPs are mainly 
protons, electrons and a particles with energy typically ranging from 10 to several hundreds of 
MeV. OCRs generally originate from outside the Solar System, e.g. in supernova remnants, and 
their composition consists mainly of protons, ~ 7-10% helium and ~ 1 % heavier elements. Due 
to the scattering effect of the magnetic fields in interstellar space, charged OCRs are subject to 
continuous deflection and the observed spectra are mostly isotropic. The GCR flux in t he Solar 
System is inversely modulated by the variations of the solar activity (e.g.. |Parkerlll958h . In the 
long term, during solar activity maximum the increased solar and heliospheric magne tic fields are 


more efficient at hindering low -energy OCRs from entering the inner heliosphere (e.g.,|HeberetaL 


20071: IWibberenz et al.l 1200211 than at sol^ act i yity minimum when the interplanetary magnetic 
field strength are reduced (iGoelzer et al.l 12013l: ISmith et al.l 12013l: IConnick et al.l l20 1 lb . Conse¬ 
quently, the GCR pop ulation is most intense during solar minimum (e.g., iMewaldt et al.l 1201 Ol: 
Schwadron et al.l 1201211 . In the short term, the GCR spectrum can also be altered indirectly by 
solar events such as coronal mass ejections (CMEs) where the enhanced interplanetary magnetic 
field can s weep away a fraction of OCRs causing reducti ons in GCR doses in the form of Eorbush 
decreases ( Eorbush 1938 : Schwadron et al. 2012 . boMhIl . 


The Mars Science Eaboratory (MSE) spacecraft (IGrotzinger et al.ll2012h . carrying the Curios¬ 
ity rover, was launched on November 26, 2 011 and the rover lan ded on Mars on August 6, 2012. 
The Radiation Assessment Detector (RAD, iHassler et al.ll2012h onboard is an ene rgetic particle 


detec tor and carried out radiation measurements during its cruise from Earth to Mars (IZeitlin et al. 


2013h and now continues to do so on the surface of the planet (IHassler et al.l 1201411 . The solar 
modulation of the dose rate measured by RAD during the MSE cruise phase has been discussed by 
Guo et al. ( 2015b l. 


Typically, CME and SEP events are more common during solar maximu m, but the latest 


solar maximum has been a very weak one compared t o space-age averages (e.g.. lSchwadron et al. 


McComas et al. 2008 


201 ll: iMcComas et al.ll2013l:lKomitov & Kaftanll2013h. This may have been caused by the reduced 


Schwadron et al.ll2014ah in the deep cycle 23 


solar wind pressure (e.g., 

24 minimum which has allowed the ter mination shock to m ove closer to the Sun and led to a 
weakened modulation of the heliosheath ( Scherer et aklboi 1 1. There were only three solar particle 
events detected by RAD over its first Martian yearQ on the surface of Mars. They can be seen in 
the middle panel of EigurefUas three peaks of the measured dose rate (the dose rate measurement 
will be explained in more detail in Section ^ at sols 242, 420, and 504 respectively. The dose 


' 1 Martian year a; 668 sols; 1 sol = 1 Mars day a; 1.03 earth day 
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rates during these SEP events were several times higher than the quiet-time dose rate. For a eloser 
inspeetion of the variations of the GCR-driven dose rate during the solar quiet periods, we omit 
the peak values of the SEPs in this figure. A zoomed-out figure eontaining the surface dose rate 
measurements for t he first 300 sols, as well as the peak doses of the first SEP can be found in 
Hassler et al. ( 2014 ). 


Th e bottom panel of Figure |T] shows the Mars surface atmospheric pressure as measured by 
REMS (iGomez-Elvira et al.ll2012l) . while the middle one shows the dose rate measured by RAD. 
Both panels show high time resolution data as a shaded band and one-sol averages as solid dots. 
The variations in the shaded bands are real short-term (diurnal) oscillations in dose rate and are 
anti-correlated with corresponding diurnal pressure variations . A zoomed-in figur e containing one- 
sol variations of the dose rate and the pressure can be found in iRafkin et al .1(120141) . However, some 
of the long-term evolutions, e.g., the drops in both dose rates and pressure in the time period pro¬ 
ceeding sol 200 to those around and after sol 300, can not be well explained by the anti-correlation 
between them. The aim of this paper is to investigate the relative influences of atmospheric pres¬ 
sure (bottom panel) and heliospheric modulation (top panel) on the measured dose rate (middle 
panel). 

During solar-quiet periods, several factors on different time scales may affect the variation of 
the GCR-induced dose rates measured by RAD on the surface of Mars: [a] the short-term diurnal 
variations of the Martian column mass (me asured as surface p ressure in a hydrostatic atmosphere) 
at Gale crater caused by daily thermal tides ( Rafkin et al.ll2014 ). [b] the long-term seasonal changes 
of the atmospheric pressure shown in the bottom panel of Figure [H and [c] the modulation of 
the primary GCR flux by the solar magnetic field which correlates with the solar activity and 
the rotation of the helio sphere. A commonly us e d parameter of helio spheric modulation is the 
modulation potential O (iGleeson & Axfordlll968l: lUsoskin et al.ll2005h which corresponds to the 
mean electric potential that quantifies the energy loss of a cosmic ray particle experiences inside the 
heliosphere and is often used to parametrize the modulation of the GCR spectrum. The modulation 
potential allows specification of distributions across a range of GCR particle species with different 
nucleons A and charge state Z and has been often used in determining GCR spectra and flux 
based on analytic models (e.g.. lONeilll 120061) . In other words, the GCR-driven primary flux in 
Equation [H is a function of O, i.e., Fo,(0) and the yield matrix, M,y, is a function of atmospheric 
pressure, P\ the resulting particle spectra, Fj(E), along with the dose rate which is a good measure 
of the radiation environment, are consequently a function of O as well as pressure, i.e. Fj{E) = 
Fj{(t>,P,E). This study aims to derive an empirical expression for the dose rate as a function 
of O and pressure based on observational data and thus make it possible to predict the radiation 
environment on the surface of Mars under different solar modulations and pressure variations. 
From the neutron monitor count rates {CRnm, in the unit of counts per minute) recorded by the 
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Fig. 1.— Top: Solar modulation potential O at Earth derived from Oulu neutron monitor count 
rate is shown in black and the per-sol-averaged O corrected to Mars’ location is plotted as red dots. 
Middle: Dose rates recorded by RAD in the silicon detector B (gray curve) and plastic scintillator 
E (tissue-equivalent, black curve) with their per-sol-averaged values marked as magenta- and red- 
dots respectively. Bottom: The pressure data from REMS (gray curve) and the per-sol-averaged 
values (red dots). The overall average pressure of the first 700 sols is about 840 Pascal and is 
marked as a thick-green line. Note that we have given the unit of time in both sol (i.e., time since 
the landing of MSE) at the bottom of the figure and day of year/year format at the top of the figure. 
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Oulu neutron monitor [§, the modulatio n potential O (in the unit of MV) at Earth ean be estimated 
( Usoskin et al. 2002 : Guo et al]l2015b l and its results are plotted in black in the top panel of Figure 
[U However, the potential O at Earth and at Mars may differ due to the longitudinal difference of 


the modulation across the Parker spirals resulted from three dimensional drifts (e.g.. ljokipii et al. 


20041: IPotgieter & Le Rouxlll992h as well as the small radial gradient between 1.0 and 1.5 AU 


Schwadron et al. 

2010; 

Gieseler et al. 

2008) 


analytic function given by 


20081). Th is radial effect can be corrected following the 


Schwadron et al. 


(120101) and the resulting modulation O at Mars (per- 


sol-average) is plotted in red dots. It is generally smaller than O at Earth with variant differences 
(shown as red-shaded areas) through time due to varying distances of the planets to the Sun. 


Transient effects that are localized to narrow ranges of heliospheric longitude, such as narrow 
CMEs, can also perturb the GCR fluxes differently at Earth and Mars. In order to reduce the 
spatial longitudinal discrepancy of O, we use binning techniques in our current study as presented 
in Section 14.2.21 


2. RAD Measurements 


RAD measures both the charged as well as the neutral radiati on environment on Mars (IHassler et al. 


2012I). It uses the dE/dx vs total E or multiple dE/dx techniques dEhresmann et aklboidl: Iguo et aP 


2015ah to identify charged particles in a telescope stack of three silicon semiconductor detectors. 


A, B, and C, followed by a high-density Csl scintillator calorimeter, D. The Csl crystal together 
with a plastic (BC-432m) scintillator (namely detec tor E) are enclosed in a highly efficient antico¬ 
incidence in order to measure the neutral radiation (IKohler et al.ll2014h . In addition, the dose rate 
is measured in the sili con detector B as we ll as in the plastic detector E. A detailed overview of the 
instrument is given in 


Hassler et al.l ( 2012h . 


In this work we will determine the influence of heliospheric modulation and atmospheric pres¬ 
sure on dose rate in silicon (detector B) and in plastic (detector E). The dose rate is defined as the 
energy deposited by radiation per unit mass and time and is measured in Gy/day (J/kg/day). During 
quiet times, the dose rate on the surface of Mars is - apart from a very small natural background 
- determined by the GCR and its interaction with the atmosphere and soil. It can be described by 
the following equation: 


E>(0,E) 


^max^ 


/l/(E, 6)E;(0, P, E)dEde/m, 


( 2 ) 


^The Oulu count rate data have been obtained from http;//cosmicrays.oulu.fi/ and the pressure effect has been 
corrected. 
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where E is the energy of a particle with type j, F/O, P, E) (in the unit of counts/MeV/sec/cm^) 
is the surface particle spectrum (equivalent to Ej{z, t, E) in Equation [T]) which is modulated by the 
heliospheric potential, O, and atmospheric pressure, P. m (kg) is the mass of the detector and 
6 is the energy deposited by the particle i n the detecto r. This energy deposit can be estimated 
using either a simple Bethe-Bloch A nsatz (lBethelll932l) or with more sophisticated models such 
as GEANT4 ( Agostinelli et al. 2003 ) and HZETRN ( Wilson et al. 1995 ) and is included as a yield 
matrix, /!/(£', e), in the above equation. 6„,;„ and ef„ax are the minimum and maximum energy over 
which the detector is sensitive and D is the corresponding dose rate integrated over the entire 
detection area (area) and all the detected particle species. Dose rate is in the unit of MeV/kg/sec 
and can be transfered to qGy/day. Correspondingly D depends on both heliospheric potential, O, 
and on atmospheric pressure, P. 

RAD measures radiation doses induced by both charged and neutral energet ic particles in two 
detectors: the silicon detector B and the plastic scintillator E (IHassler et al.ll2012l) . RAD is directly 


mounted on the ’shoulder’ of the rover body (IHassler et al. I 2 OI 2 I) and the shielding of the rover 
from above can be ignored. It is indeed shielded by the rover from below. However since the 
upward flux (of albedo particles) is much smaller than the downward flux[^ the measured dose can 
be roughly assumed equivalent to the radiation dose of the Martian surface. 


RAD operates on adjustable ’’observation” cycles, with typical durations of 32 minutes early 
in the surface mission and 16 minutes later in the mission to date. These cadences typically yield 44 
or 88 measurements per day or sol for dose rates or charged/neutral particle count rates. Detector E 
has a composition similar to that of human tissue and is also more sensitive to neutrons than silicon 
detectors. Eor a given incident flux, the dose rate in detector B is generally less than the dose rate 
in E because of the comparatively larger ionization potential of silicon as shown in Eiguredl 


3. Martian Atmospheric Pressure 


The Martian atmosphere is roughly 1% as thick as that of the Eart h’s. The dominant compo¬ 
sition of the Martian atmosphere is about 95% CO 2 (lOwen et al.lll977h . of which 25% condenses 
seasonally onto the winter pole. Nevertheless, the seasonal variation of the composition has a 
very little effect on the surface radiation field compared to the changes of the column mass (e.g., 
Rafkin et al .1120141) . The oscillations in atmospheric column mass drive variation of energetic par¬ 
ticle radiation at the surface with both diurnal and seasonal periods. 


^More detailed analysis shows that the upward-downward ratio of the proton flux is only about 13% (Appel et al., 
paper in preparation). 
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3.1. Diurnal Pressure Variations 


With a constant value of gravitational acceleration g, the surface pressure is an exact measure 
of the column mass given a hydrost atic atmosphere. A n increase in column mass corresponds to 
an increase of surface pressure (e.g-. lRafkin et alJl2ni4t) . 


The Martian atmosphere exhibits a strong thermal tide excited by direct solar heating of the 
atmosphere on the dayside and strong infrared cooling on the nightside. Heating causes an inflation 
of the atmosphere with a simultaneous drop in surface pressure. In Gale Crater, the thermal tide 
produces a diurnal variation of column mass of about ±5% relative to the rnedian , as measured by 
the Rover Environmental Monitoring Station (REMS) (IGomez-Elvira et al.ll2012h . The magnitude 
of the diurnal pressure cycle at Gale Crater is substantially greater than previous surface measure¬ 
ments. This is likely due to the topography of the crater enviro nment, which yields hydrostatic 
adjustment flows that amplify the daily tides (IHaberle et al.ll2014h . 


This daily pressure oscillation can be seen in the bottom panel of EigurefH where the pressure 
(shown as gray lines) expands away from the daily averaged pressure (shown as red dots) within 
the range of « + 50 Pascals. 


3.2. Seasonal Pressure Changes 


The seasonal Martian atmospheri c pressure variation is contro lled by a complex balance be¬ 
tween the cold and warm poles (e.g., iTiflmanll 19881: IZureklll988h . As on Earth, when the south 
pole is in total darkness, the north pole is experiencing continuous sunlight; one might expect that 
the global pressure should stay roughly constant over the year, as CO 2 vaporized at one pole would 
freeze at the other. However, the high eccentricity of Mars’ orbit causes the insolation to be signif¬ 
icantly different between poles. Mars is farther from the Sun during northern summer; the summer 
in the southern hemisphere is much warmer than summer in the northern hemisphere. As a result, 
the north and south poles have different impacts on the atmospheric pressure changes through con¬ 
densation of CO 2 to the polar region in winter and the recession of CO 2 polar cap during spring 
and summer. The seasonal CO 2 condensation cycle results in the seasonal pressure variation: as 
more CO 2 evaporates into the atmosphere in the summer, the measured surface pressure increases 
as shown in the bottom panel of EigureH] 


Curiosity landed at the time when the northern hemisphere was in late summer and the global 
pressure was near its minimum since the southern CO 2 ice cap had nearly reached its maximal 
extent during southern hemisphere winter. As shown in Eigure [B the atmospheric pressure then 
began to increase as the southern polar cap started shrinking during the southern spring (northern 
autumn). It then reached a peak (~ sol 175) during early northern hemisphere winter, when the 
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southern cap was near its minimum level and before the northern eap had grown to its maximum 
size. A small minimum of the pressure is present during the late northern winter when the northern 
eap reaehes its maximum extent (~ sol 310). 

The diurnal pressure oseillation (as described in Section IXTI) superimposed on this long term 
seasonal pressure change ean be averaged out by calculating the per-sol-averaged pressure as 
shown by the red dots in Figure \T\ The peak to peak seasonal pressure difference reaehes about 
25% of the average pressure. This seasonal pressure variation, like the diurnal variation, eauses 
eolumn density changes of the atmosphere that affeet the particle fluxes measured by RAD at a 
seasonal period. 


4. Empirically Modeling the Dose Rate as a Function of Pressure and Solar Potential 


As discussed in Seetion [TJ we aim to find out how RAD measured surface radiation environ¬ 
ment depends on atmospheric pressure P and solar modulation O. Both influenees are blended and 
embedded in the long-term variations of the measured GCR dose rate as seen in Figure [B While 
some previous studies have attempted to model the effeets of atmospherie pressure on the Martian 
radiation environment and dose rate0, this work has the advantage of using aetual in-situ measured 
data for a quantitative and empirieal study. 


Rafkin et al.l (l2014t) have analyzed the diurnal pressure effeet on the dose rate measurement by 


successfully isolating the diurnal variations in the RAD measurements from the longer-term influ¬ 
ences that include seasonal atmospheric shielding and variability of the heliosphere. This method 
will be explained in Section 14.1.11 The authors presented a robust fit with a linear correlation be¬ 
tween pressure ehanges and dose rate measurements @ providing the reaso nable assurnption th at the 
solar modulation O varies at a time seale longer than the diurnal period. IGuo et al.l (l2015bh have 
studied the solar modulation of the GCR dose rate measured by RAD during the MSL eruise phase 
when the pres sure-variation effeet was not present. Two separate empirieal models were employed 
to deseribe the anti-eorrelation between heliospherie modulation potential O and measured dose 
rate. Both a simple linear funetion and a non-linear regression model eould equally well represent 
this anti-eorrelation given that the shielding of the spaeeeraft (’pressure’) was constant. Assuming 
that O and pressure P are two independent parameters infiueneing the surface dose rate (while the 
viability of this assumption will be discussed in Section [ 5 ]), the pressure effect and the solar modu- 


lEhresmann et all j201 ih calculated the surface radiation exposure for much higher atmospheric pressures which 
might have been present during the Noachian epoch. 

^This linear correlation has been obtained at the scale of pressure values measured at Gale Crater and should not 
be simply extrapolated over a much wider range of pressure, e.g., to the top ot the Martian atmosphere. 
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lation can be linearly combined. This results into the following two models of GCR-induced dose 
rate variations on the surface of Mars. 


1. Both the pressure effect and the solar modulation drive the variation of dose rate linearly and 
independently, written as: 


p) - Dqi + kP +^o, 


( 3 ) 


where k (in the unit of p,Gy/day/Pa) is the linear correlation factor between pressure and dose 
rate variations and can be fitted when O is constant; yS (in the unit of [iGy/day/MV) is the 
linear correlation factor between solar potential and dose rate changes and can also be fitted 
when pressure is approximately stable; Dqi (in the unit of p,Gy/day) is some relevant dose 
rate when both pressure and O are at certain typical levels, e.g. Pq and Oq. 


2. The so lar modulation effect can also be described by a non-linear e mpirical model (IGuo et al. 


2015bh as often used in the analysis of neutron count rates (e.g.. lUsoskin et al .11201 It) . The 


combination of the pressure and O changes results in: 


D{Q>,P) = Dq2+kP + -^ 

(D + a2 


( 4 ) 


with the requirement ofa 2 > 0 to assure a positive denominator. ai (in the unit of MV • Gy/day) 
and a 2 (in the unit of MV) can be fitted when pressure is constant. Dq 2 (in the unit of 
qGy/day) is some relevant dose rate when both pressure and O are at certain typical levels. 


4.1. Pressure Effect 


The day and night variations of the dose rate have bee n observed for the first time on the 
surface of Mars by RAD, as presented bv iRafkin et al.l (120141) where the anti-correlation between 
pressure and dose rate changes has been quantitatively investigated. Due to the day and night 
oscillations of the atmospheric column mass, the characteristics of the particle radiation at the 
surface also vary diumally. The middle panel in Figure [H shows the dose rate measured by RAD 
on the surface of Mars during the first 700 sols (06/Aug/2012 - 25/Jul/2014). The dose rate varies at 
a diurnal level seen as oscillations of the black curve for the plastic detector, E, and the gray curve 
for the silicon detector, B. The variation in the E dose rate over a diurnal cycle averages about 15 
qGy/day (or ~5 %) peak to peak, out of around 220 qGy/day. The oscillation in the silicon detector 
is more pronounced because it includes both the diurnal variation and larger statistical fluctuations 
due to its smaller geometric factor. The per-sol-averaged values of the dose rates are calculated and 
plotted over the curves as magenta-dots (for B measurements) and red-dots (for E measurements). 
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Three SEP events occurred on sol 243,421 and 504 and they are excluded from the following study 
of the pressure effect since the particle fluxes and energies during SEPs are substantially different 
from those of OCRs and the pressure response of dose rate can be heavily modified. Quantitative 
analysis has shown a clear inverse relation between the variations in the atmospheric pressure and 
the RAD dose rates with a correlation coefficient for linear regression of 96%. 


4.1.1. Fitting k using Hourly Perturbation 


We use the method described in lRafkin et al .1(120141) to produce the average diurnal perturba¬ 
tions of the data; this approach aims at isolating the diurnal variations in the RAD measurements 
from the longer-term influences that include seasonal atmospheric shielding and variability of the 
heliosphere, as well as solar event. 


Generally, the dose rates measured by RAD are distributed uniformly in time while REMS’ 
pressure data are recorded at 1 Hz over 5 minute periods at certain periods of the sol. In order to 
correlate these two data sets with different time frames, we obtain the hourly averaged measure¬ 
ments of both pressure and dose rates and their corresponding hourly perturbations. The hourly 
averaged dose rate and pressure are D/, ^ and respectively where h represents the hour (here 
defined to be 24 hours on each sol) and s corresponds to the sols. The measurements, especially 
the pressure data, were not always uniformly taken in time and there are gaps in the data lasting 
longer than one hour. Therefore we only consider sols where there are at least 20 hours with pres¬ 
sure measurements in each hour. Then we interpolate the pressure data using a spline interpolation 
method so that data are distributed uniformly and continuously over 24 hours through that sol and 
systematic errors of the hourly average and per-sol-average could be minimized. Einally the in¬ 
terpolated pressure data can be binned into 24 bins for each sol and the hourly binned pressure 
is the corresponding hourly average Ph^s- The hourly pressure perturbation dR/, ^ is defined as the 
difference between the hourly average, Ph,s, and the total average, Ps, of the corresponding sol: 

24 

Z Ph.s 

SPh.s=Ph,s-Ps = Ph..-—^. (5) 

In this case, the hourly perturbation is isolated from the sol-to-sol changes; 6Ph,s can further be av¬ 
eraged through all sols to obtain the mean hourly perturbation 6Ph and the corresponding standard 
error in that hour. The same binning technique has also been applied to RAD dose rate measure¬ 
ment. The mean hourly perturbation of dose rate, dD/,, can be readily correlated with the hourly 
pressure perturbation, 6Ph, and their relationship follows a clear anti-correlation which can be fitted 
with a first-order polynomial function: 


6Dh = KdX 5Ph. 


( 6 ) 
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4.1.2. Pressure Ejfect on Dose Rates 


Using the binning method deseribed in Seetion l4.1.1[ we fitted the diurnal dose rate variations 
(measured by deteetor E) as a response of pressure oseillations through the first MSL Martian year 
as shown in Figure[2l The eorrelation eoeffieient for linear regression is 98.8%. The fitted propor¬ 
tionality faetor, k, is -0.13 + 0.02 p,Gy/day/Pa where the error bar is obtained through propagation 
of the uneertainties of the hourly perturbation data. The measured dose rate deereases when the 
pressure inereases sinee k is negative. This indieates that the atmosphere has a shielding effeet on 
surfaee dose rates . The proportionali ty faetor obtained here is eonsistent, within error bars, with 
what obtained by iRafkin et al.l (120141) whieh was about 0.116 tiG y/day/Pa. The slight differenee 
is eaused by the seleetion of different periods (IRafkin et al.l (120141) used the first 350 sols of data), 
and by the method used here to interpolate data into a uniform distribution in time before the bin¬ 
ning proeess. We have applied the same method to fit the pressure effeet on dose rates deposited 
in deteetor B and found the eorrelation eoeffieient to be 96.4% and the proportionality faetor to 
be -0.10 + 0.06 p,Gy/day/Pa whieh is smaller than that for the E deteetor sinee dose rates in B are 
generally smaller. 


Note that eompositional ehanges in the atmosp here have only a y anishingly small effeet on 
surfaee radiation eompared to the pressure variations ( Rafkin et al.ll2014 ). Similarly, the oeeasional 
presenee of dust in the atmosphere during dust storms is negligiblql 


4.2. Solar Modulation 

The pressure variations and the modulation of the primary GCR radiation outside the atmo¬ 
sphere are the two faetors whieh determine the longterm dose rate variations on the surfaee of Mars 
as explained in Seetion In order to separate these two effeets, we subtraet the pressure effeet 
using the proportionality of the variations found above. 


4.2.1. Subtracting the Pressure Ejfect 

As we ean already isolate the pressure effeet and obtain the proportionality faetor k in Equa¬ 
tion [3] and S the effeet of solar modulation ean be investigated by first subtraeting the pressure 


® A PLANETOCOSMIC simulation has been carried out to derive the particle fluxes on the surface of Mars con¬ 
sidering the existence of dust storms. The result shows that the effect of dust storms is very small (Appel et al., paper 
in preparation). 
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S Pressure [Pa] 


Fig. 2.— Hourly perturbation of dose rate SD/, versus hourly perturbation of pressure 6P/i through 
sol 13 to 682 (approximately one Martian year) as shown in blue. The error bars stand for the 
standard deviation of the averaged hourly perturbation. The fitted anti-correlation is shown as a 
red line with a slope of Kd being -0.1306 + 0.0176 [iGy/day/Pa. 
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Fig. 3.— RAD dose rates (pGy/day, left axis) and solar modulation potential O at Mars’ radial 
distanee (MV, think-blaek-dotted-line, right axis) through sol 13 to 682 (approximately one Mar¬ 
tian year). RAD dose rates from silieon deteetor B is in blue and from plastie deteetor E is in red. 
Thiek-green and thiek-eyan lines are the resulting dose rates {D' in Equation|7]) assuming pressure 
is eonstant at Pq. Note that we have given the unit of time in both sol (i.e., time since the landing 
of MSL) at the bottom of the figure and day of year/year format at the top of the figure. 
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efFect from dose rate and then eorrelating the ’eonstant-pressure’ dose rate with solar modulation 
potential O. The method can be described by restructuring, for instance Equation[3]as following: 

D'(0) = D'i+;S0, (7) 

where D' = D - k{P - Pq) and = Dqi + kPq. Employing all the per-sol-averaged dose rate 
data collected over one Martian year (sol 13 to 682), we can subtract the seasonal pressure effect 
assuming the surface pressure is constant at a particular value, i.e., Pq = 840 Pascals, which is the 
averaged pressure found by REMS during this time period as shown in Eigure [B By calculating 
the difference of the per-sol-averaged pressure P and Pq and using k found in Section I4.1.2[ we 
estimate the pressured-induced dose rate to be k{P-Pq) which is then subtracted from the long-term 
dose rate measurement with the ’constant-pressure-as-Po’ dose rate remaining, namely D' being 
solely a function of solar modulation potential O. Eigure [3] shows the RAD dose rates measured by 
both silicon and plastic detectors before (thin lines) and after (thick lines) the pressure correction. 
Also plotted is the solar modulation O which already shows an anti-correlation with the constant- 
pressure dose rates. 


4.2.2. Correlation of dose rates and O 


The modulation potential, O, is derived from the Oulu neu tron monitor on Earth (a t lAU) and 
corrected to the radial distance of Mars (~ 1.5 AU) following ISchwadron et al.l (l2010h . However 
Mars and Earth are not always magnetically well connected and cross-field diffusion and drift can 
be extremely important. In other words, the modulation process is fundamentally 3-dimensional 
and O can not directly represent modulations at Mar^. In order to average out the cross-field 
discrepancy of the modulation between Earth and Mars and to smooth out the rotation of the 
heliospheric magnetic fields, we bin the per-sol-averaged data as shown in Eigure [3] into 26-sol- 
averaged bins and correlate the binned O values and constant-pressure dose rates. The correlation 
coefficients are -0.66 and -0.56 for detectors B and E respectively, clearly indicating a negative 
correlation between solar modulation and GCR induced dose rates. We applied regression fittings 
of both the linear and non-linear models described by Equations[3]and|4]to data from both detectors 
as shown in Eigure IH 

In order to reliably propagate the uncertainties due to bot h measureme nts and binning pro¬ 
cesses, we have carried out a bootstrap Monte Carlo simulation (lEfronlll981r) . Simulated data sets 
are generated using the uncertainty range of the binned data which are then fitted to the models; 


^During the cruise phase MSL was mostly magnetically connected with the Earth (IPosner et alj|2013l : IGuo et al. 


2015bh and directly correlating <I) measured at Earth and the RAD dose rate was sensible. 
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Fig. 4.— 26-sol binned data and fittings processed through bootstrap Monte Carlo simulations of 
the variations of the RAD dose rate (with pressure assumed to be constant Pq = 840 Pa) and the 
solar modulation O on the surface of Mars through sol 13 to 681.The red/magenta line and area 
represent the linear fits with standard errors (Eq. |7]) of constant-pressure dose rate in detector E/B 
versus solar modulation potential. The green/cyan line and area show the results of the non-linear 
fits (Eq.lH). 
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dose rate results at a wider range of O are estimated simultaneously at eaeh fit; 500 simulated fits 
were proeessed for eaeh model; for every O value ranging from 250 to 1200 MV, 500 dose rate 
values were generated and their mean is taken to be the ’predieted’ dose rate while their standard 
deviation as the uneertainty. The fitting results of the two separate models applied to both deteetors 
are presented below. 

• The fitted parameters jS and for the linear model are obtained as the mean values of the 
500 fitted parameters and their uncertainties are the standard deviation of the 500 Monte 
Carlo fits. For the silicon detector, B, we obtained the parameters of the linear model as 
Dqj = 224.10+16.12 p,Gy/day andyS = -0.12+0.03 p,Gy/day/MV. Forplastic detector E, the 
fitted parameters are = 267.17 + 16.29 p,Gy/day andyS = -0.11 +0.03 p.Gy/day/MV. Note 
that the absolute values of are not essential to our study and they depend on the choice 
of normalized pressure, Pq. The parameter yS however, as shown in Equation |71 directly 
represents the linear dependence of dose rate changes on the solar modulation potential, O. 
The results obtained from measurements of the two different detectors are consistent with 
each other within error bars. 

• The non-linear model has three fitting parameters: a\, a 2 and 0^2 = ^02 + /ePo- Their values 
and error bars are also obtained using the same Monte Carlo method. Eor silicon detector, 
B, the fitted parameters are Z)q 2 = 89.4 + 14.0 MV, ai = (3.7 + 0.8) x 10^ MV [iGy/day and 
ai = (9.9 + 0.1) X 10“^ MV. Eor plastic detector, E, the results are 0^2 = 149.6 + 15.3 MV, 
ai = (3.3 + 0.8) X 10"^ MV jrGy/day and a 2 = (1.3 + 0.1) x 10“^ MV. The small absolute 
values of a 2 indicates that it could be ignored in the model for the range of data in the current 
measurement. 

The ’predicted’ dose rates (with constant surface pressure Pq) at given O values from 250 MV 
to 1200 MV were estimated for both models and are plotted in EigurelH The uncertainty of the 
predicted dose rate increases when the extrapolation is further away from the actual measurements. 
The linear model often predicts a smaller dose rate than the nonlinear model and this difference 
is bigger for small solar potentials, i.e., during solar minimum. Eor instance, at O = 250 MV the 
discrepancy between the two models is as large as 80 qGy/day for the silicon detector, B. Because 
the predictions of the two models differ substantially at solar extreme conditions, choosing one 
or the other for predicting the radiation exposure of an astronaut would make a big difference. 
The data themselves do not rule out one or the other because the time period for which they are 
currently available does not cover a sufficiently large range of O. This underlines the importance 
of acquiring more data over an extended period of time to adequately cover the entire solar activity 
cycle. 
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4.2.3. Estimates of the Surface Radiation Environment under different O and P 


The annually av eraged modulation potential O reeonstructed from Oulu neutron monitor data 
(lUsoskin et alJl201 ih from 1937 until 2014 ean be used to calculate the corresponding ’expected’ 
dose rate predicted by our models. We can also use the sunspot number predicted by the U.S. Dept, 
of Commerce, NOAA, Space Weather Prediction Center (SWPC) for the coming years (2015 to 
2019) to estimate the corresponding solar potential following a correlation study of monthly O and 
sunspot numbers (IGuo et al.ll2015bh . The annual modulation potential, O, extrapolated to 1.5 AU 
radial distance from the Sun through year 1937 until 2019 is shown in the top panel of Figure[5l It 
shows a clear 11-year cycle and varies from less than 250 MV to more than 1200 MV. 

For evaluating the space radiation environment, dose equivalent is often derived and can be 
assumed to be proportional to the risk of lifetime cancer induction via population studies (ICRP60). 
Its approximated value, in Sieverts (Sv), is taken to be ,< Q > x D , where D is measured tissue- 
equivalent dose and < Q > is the average quality factor which is a conventional parameter for 
radiation risk estimation. A < Q > value o f 3.05 + 0.26 was fo und from RAD’s measurements on 
the surface of Mars over the first 350 sols (IHassler et al.ll2014t) . Multiplying the tissue-equivalent 
dose rates (directly measured by the plastic detector E or modeled by e.g. HZETRN model) with 
the average quality factor yields the dose equivalent rate of OCR fluxes. Eor d ose rates measured 
by detector B, a silicon to water conversion factor of 1.38 has to be applied first ( Zeitlin et al.ll2013 ) 
which approximately relates energy loss per unit of path length {dE /dx) in silicon to Einear Energy 
Transfer (EET) in water. RAD measured dose equivalent is thus most comparable to high-water 
content skin dose equivalent. The body effective dose can be further derived as the weighted sum 
of different organ dose equivalent (skin, eye, bone, brain, heart, etc.) and the weighting factor 
for each organ can be foun d in the National Council on Radiation Protection and Measurements 
( Einton & Mettler Jr 2003 ). 


The extrapolated dose rates at different modulation potential values have been estimated via 
two different models and two different detectors as shown in EigurelH These four sets of modeled 
dose rates can be transfered into dose equivalent rate and the results mostly agree with each other 
within error bars §. Because we can not currently determine which of the two models is the better 
approximation when the solar potential is outside the measured range, we let all the four modeled 
values serve as possible results and take their mean as predictions and the propagated errors as 
uncertainties. The thus final modeled dose equivalent rates and their standard deviations are shown 
in Eigure[5l 


^Figure 4 in 


Guo et at 


(l2015bh has shown the modeled dose equivalent rate during the cruise phase estimated by 


both linear and non-linear models based on dose rates from both detectors. The results are consistent within error bars 
with exceptions during extreme solar conditions for which the non-linear model predicts higher dose rates. 
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Fig. 5.— Top Panel: Annual average of reeonstructed modulation potential O at the radial distanee 
of Mars (gray dotted line) sinee 1937 until 2014 (on the left side of the vertical dashed line) 
as well as the predicted values from year 2015 to 2019 (on the right side of the vertical dashed 
line). Bottom Panel: Modeled annual dose equivalent rate following the evolution of O assuming 
different seasonal surface pressures: blue for po = 700 Pa, green for po = 970 Pa and red for 
Po = 840 Pa. 
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Note that the modeled dose equivalent rate has been derived based on the ’eonstant-pressure- 
at-840 Pa’ dose rate (see Seetion l4.2.1l) . However, the seasonal surfaee pressure at Gale erater, as 
shown in Figure [B expands between 700 and nearly 1000 Paseals whieh would result in roughly 
0.1 mSv/day of dose equivalent rate ehanges§ and this is a eonsiderable portion ~ 16% of annual 
averaged surfaee dose equivalent rate ~ 0.7 mSV/day. In order to show this seasonal pressure 
effeet, we also estimated the dose equivalent rate (at various solar potentials) when pressure is 700 
and 970 Pa respeetively as shown Figured The dose equivalent rates are inversely related to the 
surfaee pressure, although the seasonal pressure influence is much smaller than the longer term 
effects driven by solar modulation. 

The estimated surface dose equivalent rate ranges from about 0.35 mSv/day to about 1.15 
mSv/day and has a clear anti-correlation with O as expected from both models. Stronger solar 
modulation leads to a decrease of dose equivalent rate and at O > 1000 MV the dose equivalent 
rate can be as low as ~ 0.35 mSv/day within the uncertainties, considerably smaller than the cur¬ 
rent averaged measurement. At solar maximum and minimum when the modulation potential is 
further away from the measured range, the uncertainty of the estimations increases due to the large 
discrepancy between the models. Future measurements over solar minimum periods are essential 
for improving the predictions at low modulation potentials. 


Due to the shielding of the atmosphere, the current su rface dose equivale nt rate is only 
about 40% of the RAD cruise measurement ~ 1.8 mSv/day dZeitlin et al.l 120131) and that from 
the Cosmic Ray Telesco pe for the Effects of Rad iation (CRaTER) on the Lunar Reconnaissance 
Orbiter ~ 1.6 mSv/day (ISchwadron et al.ll2014bh . Consequently, our estimations of the surface 
dose equivalent rate over solar minimum and maximum periods are muc h smaller than the pre¬ 
dictions b ased on deep space mea surements given by RAD (cruise phase. IGuo et al.ll2015bh and 
CRaTER (ISchwadron et al.ll2014bt) . The predictions of dose equivalent rates behind these space¬ 
craft shielding conditions are between about 1 mSv/day (solar maximum) and 5 mSv/day (solar 
minimum). 


5. Discussion and Conclusions 


We have presented the dose rate data collected by MSL/RAD on the surface of Mars and 
analyzed its short-term and long-term vari ations which ^e dr iven by the atmospheric pressure 
changes and solar modulation. Eollowing iRafkin et al.l (120141) . we first analyzed the dose rate 
dependence on diurnal pressure oscillations over the first MSL Martian year. This pressure-driven 


®250 Pascals of pressure change would lead to 32.5 pGy/day of dose rate difference which is then transfered to 
dose equivalent rate via the quality factor. 




















-21 - 


effect on dose rate changes in the longterm is then subtracted to recover the constant-pressure 
dose rate which is then correlated with solar modulation potential O. A clear anti-correlation, 
with a correlation coefficient of about -0.6 between O and the recovered dose rate, suggests that, as 
expected for higher solar activity, GCR particles are more attenuated and the dose rate is decreased. 


We carried out a quantitative study of this anti-correlation and fitted two models to the mea¬ 
sured data using a bootstrap Monte-Carlo method to estimate the uncertainties. The predictions 
of the two models for solar activity minimum differ substantially and th e data are insufficient to 
decide which of the two models should be used. ISchwadron et al.l (l2014b[) have estimated the deep 
space dose rate at different modulation potential derived from HZETRN model and the dose rate 
is indeed non-linearly dependent on O. ffowever the shape (parameters) of this analytically de¬ 
rived model is different from that of our empirically modeled function. This highlights the need 
for extended measurements to cover solar activity minimum and an entire activity cycle and these 
observations can be used to constrain the analytic models. 

The extrapolated dose equivalent rate at various modulation potentials and different surface 
pressures are shown in Figure [51 The predicted dose equivalent rate during solar maximum years 
(e.g year 1991) when O ~ 1100 MV was found to be as low as ~ 0.35 mSv/day, which is consider¬ 
ably lower than the current surface measurement ~ 0.7 mSv/day since the current solar maximum 
is atypically quiet. The modeled dose equivalent rate under solar minimum conditions can be as 
high as 1.15 mSv/day within the uncertainties. The seasonal pressure changes may affect the esti¬ 
mated dose equivalent rate at a level of about 0.1 mSv/day. Although this is less than the long-term 
solar modulation effect, it should not be ig nored. Based on the solar modulation potential predicted 
for the next five years (IGuo et al.ll2015bh . we estimate a trend of increasing dose equivalent rate 
(between 0.56 and 0.84 mSv/day) from 2015 until 2020. 

The correlation between hel iospheric potential and RAD-measured dose rates during the 
cruise phase was investigated by IGuo et al.l (1201 5bt) with the data modeled by both linear and 
non-linear functions. The linear dependence of dose rate on O was -0.39 + 0.07 p,Gy/day/MV 
for the silicon detector B and -0.44 + 0.04 p,Gy/day/MV for the plastic detector E. In comparison, 
yS derived here for the surface case are much smaller: -0.12 + 0.03 and -0.11 + 0.03 p,Gy/day/MV 
for B and E respectively. This is because [a] the magnetic connection between solar modulation 
potential O measured at Earth and dose rates evaluated during the cruise was very good, allowing 
the direct and thus stronger correlation of the daily values; and more importantly: [b] low energy 
particles which are more affected by solar modulation make a bigger contribution to the GCR doses 
detected during the cruise phase than to the doses measured at the surface of Mars. The process of 
a GCR spectrum penetrating thro ugh the Martian atmosp here can be simulated using for instance 
the PEANETOCOSMICS toolkit ( Desorgher et al. 2006 ). Given typical Martian atmospheric con¬ 
ditions, we have found that protons with energies less than 170 MeV do not reach the bottom of 
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Gale Crater. The penetration energy is speeies dependent and inereases with inereasing ion charge. 
Instead, the spacecraft shielding was highly non-uniform while nearly 50% of inco ming particle 


trajec tories within the field of view of RAD from above were only lightly shielded dZeitlin et al. 


20131), thus allo wing more low-energy particles to contribute to the dose rate. We will carry out 


more GEANT4 dAgostinelli et al.ll2003l) and PLANETOCOSMICS simulations as well HZETRN 
modelling in the future to study the interactions of GCR spectra with different shielding conditions 
(spacecraft and the atmosphere) in order to derive modeled function of the dose rate dependence 
on pressure and solar potential which can be compared with our empirical functions. 


Both models have assumed the independent effect of pressure and O on dose rates. However, 
this may be modified when pressure and O change over wider ranges than have been observed to 
date: 


• a much thinner atmosphere will allow more lower-energy particles to reach the surface which 
experience stronger modulation (e.g., bigger |/3| in the linear model). Therefore, for signifi¬ 
cant pressure changes, (3 could be a function of P , i.e., /3(P); 


much stronger solar modulation (bigger O) would lead to a larger fraction of high-energy 
particles in the GCR flux and these energetic particles are less affected by the atmosphere 
(smaller |a:|); when pressure is much higher and the su rface atmospheric depth is approach¬ 
ing the Pfotzer maximum (e.g.. lRichter & Raschll2008h . most primary particles are shielded 
while more secondary particles are generated and this may result in a decreased shielding 
effect; therefore the dependence of dose rate on pressure may be modified as P and O change 
substantially, i.e., k = k{P, O); 


• both the linear and non-linear models are empirical and derived from measurements; despite 
the robustness of the fitting of the actual data, the extrapolation is highly uncertain and a 
complete model requires measurements over the full range of solar conditions. 


To verify these non-linear and second-order effects of our fitted parameters in the empirical func¬ 
tion, Monte Carlo simulations as well as analytic HZETRN modeling will be carried out for con¬ 
straining the parameters and comparing the predictions. 


The quality factor used to derive equivalent dose from dose may also be sensitive to substantial 
changes in solar modulation because the estimation of t he average quality factor depends on the 
spectra of particles depositing energies in the detector (ISchwadron et al.ll2014bl) . Hitherto, the 
measured < Q > has been quite stable since [a] the RAD measurements have undergone only small 
changes of the solar modulation and [b] the modulation affects more the low-energy ions which 
are more likely to be shielded by the atmosphere and contribute little to the surface spectra. The 
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change of the quality factor under much different solar conditions needs to be investigated with 
extended measurements. 


A total Mars mission GCR dose equivalent can be e stimated based on our measurements and 
predictions of both the surface case and the cruise phase (iGuo et al.ll2015bh . The fastest round trip 
with on-orbit staging and existing propulsion technologies has been estimated to be a 1 95-dav trip 


120 d ays out, 75 days back with an extra e.g., 14 days on the surface), as described by lFolta et al. 
20121) . During solar maximum periods when O ~ 1200MV, this would result in a GCR-induced 


cruise dose equivalent of 195 + 98 mSv and surface dose equivalent of 4.9 + 2.0 mSv, which adds 
to 200 + 100 mSv during the total mission. 

Additional contributions to dose rate and dose equivalent rate by SEPs should not be ignored 
and they can differ significantly from the current measurements due to the high variability of their 
frequencies and intensities. Small, ”soft-spectrum” solar events where particle energies are modest 
will have little or no effect on the surface dose due to atmospheric shielding. For instance the total 
dose equivalent fro m the first SEP event observed on Mars (on sol 242) was only 0.025 mSv 
(IHassler et al.ll2014l) while the cruise S EP events had dose equivalent ranging from 1.2 mSv/event 
to 19.5 mSv/event dZeitlin et al.ll2013b . Future measurements of much bigger SEP events on the 
surface of Mars are crucial for understanding extreme conditions of radiation environment on Mars 
for potential manned missions. 


RAD is supported by the National Aeronautics and Space Administration (NASA, HEOMD) 
under Jet Propulsion Eaboratory (JPE) subcontract #1273039 to Southwest Research Institute and 
in Germany by DER and DER’s Space Administration grant numbers 50QM0501 and 50QM1201 
to the Christian Albrechts University, Kiel. Part of this research was carried out at JPE, Cali¬ 
fornia Institute of Technology, under a contract with NASA. The sunspot data has been obtained 
from: WDC-SIESO, Royal Observatory of Belgium, Brussels. We are grateful to the Cosmic Ray 
Station of the University of Oulu and Sodankyla Geophysical Observatory for sharing their Neu¬ 
tron Monitor count rate data. The data used in this paper are archived in the NASA Planetary 
Data Systems Planetary Plasma Interactions Node at the University of California, Eos Angeles. 
The archival volume includes the full binary raw data files, detailed descriptions of the struc¬ 
tures therein, and higher-level data products in human-readable form. The PPI node is hosted at 
http://ppi.pds.nasa.gov/, 
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